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Natural models of supersymmetry with a gravitino LSP provide distinctive signatures 
at the LHC. For a ncutralino NLSP, sparticlcs can decay to two high energy photons 
plus missing energy. We use the ATLAS diphoton search with 4.8 fb" 1 of data to place 
O j limits in both the stop-gluino and neutralino-chargino mass planes for this scenario. If the 

Q_J neutralino is heavier than 50 GcV, the lightest stop must be heavier than 580 GcV, the 

gluino heavier than 1100 GeV and charginos must be heavier than approximately 300-470 
GeV. This provides the first nontrivial constraints in natural gauge mediation models 
I . with a neutralino NLSP decaying to photons, and implies a fine tuning of at least a few 

k*" , percent in such models. 

^ '. 1 Introduction 

Supersymmetry remains the most favourable solution to the hierarchy problem. Not only 
£N| . does it elegantly cancel the quadratic divergences of an elementary Higgs boson, but it also 

achieves a precise gauge coupling unification without the need for large threshold effects. 
■ The Large Hadron Collider (LHC) has begun to systematically exclude a sizeable portion 

of the parameter space in the supersymmetric standard model. In particular, the limits 
$_i | on squark and gluino masses suggest that minimal versions of the supersymmetric standard 

model have become increasingly tuned, putting into question the original motivation for weak 
scale supersymmetry. 

However, tuning in the Higgs sector of the supersymmetric standard model depends mainly 
on the /U-term and on sparticles with large couplings to the Higgs boson, namely the third 
generation sfermions, gluinos and electroweak gauginos. If the supersymmetry breaking mech- 
anism is actually flavour dependent, third generation sfermions can be much lighter than the 
first two generations. The scenario has recently been coined "natural supersymmetry" [IHS] 
and, together with the fact that the stop production cross section is much smaller than that 
of the up and down squarks at the LHC, it allows the stringent limits on squark and gluino 
masses to be alleviated. 

A distinguishing feature of such models is whether or not the gravitino is the lightest 
supersymmetric particle (LSP). Most analyses in the literature have focussed on the case of 
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a heavy gravitino but, when the gravitino is the LSP, sparticles decay via the next-to-lightest 
supersymmetric particle (NLSP) and provide distinctive collider signals [6]. In fact there is 
a broad class of models based on new gauge forces and extra dimensions (or locality) where 
supersymmetry breaking is transmitted to the first two generation sfermions in a flavour 
blind fashion, with third generation sfermions receiving suppressed contributions via gauge 
interactions [71-119] . Some of these models can simultaneously explain the fermion mass hierar- 
chy [9~lll2| and even incorporate unification [18] but, as of yet, are not strongly constrained by 
collider searches. This motivates considering the phenomenology of natural supersymmetric 
models with a gravitino LSP, typical of models where supersymmetry breaking is mediated to 
the third generation sfermions by gauge interactions. These scenarios are collectively referred 
to as natural gauge mediation (NGM). 

We will study a particularly interesting subset of NGM models that contain a light bino, 
although our analysis also applies to other mediation mechanisms that yield a similar spec- 
trum. Natural supersymmetry is usually constrained at the LHC by searches for jets plus 
missing energy from a neutralino LSP but, in NGM with a light bino, the neutralino can 
decay to a photon and gravitino. The signal becomes two high energy photons plus missing 
energy, greatly reducing the Standard Model background, and was previously studied in [9]. 
In this Letter we recast the 4.8 fb _1 ATLAS diphoton search to find bounds on NGM models 
with a neutralino NLSP. We investigate both electroweak and coloured sparticle production, 
producing robust and otherwise model independent bounds in the stop-gluino and neutralino- 
chargino mass planes respectively. These limits provide the first nontrivial constraints on this 
class of models and allow us to quantify the degree of fine tuning required to avoid the bounds. 



2 Natural gauge mediation 

The NGM spectrum is determined by minimising the amount of fine tuning in the Higgs 
sector, retaining the features of gauge mediated supersymmetry breaking where possible. In 
particular we assume a light bino. While this is not necessary from a naturalness point of view, 
a small value for Mi is common in models with gauge mediated supersymmetry breaking. 

In a generic supersymmetric model there are three main sources of fine tuning [20]: the 
//-term 

A, = K > (!) 
m h 

one and two loop Higgs mass corrections from gauginos (with soft mass parameters M2 and 
M 3 ) 

2nm 2 & \TeV ) ' V ; 

^W^/A), (3) 



7r 3 m 2 l + tan 2 /3 V T eV 
and the one loop correction from the top Yukawa coupling (yt) 

3y 2 (m 2 Qi +ml 3 + \A t \ 2 ) tan 2 p j A ^ 



4 = 1 + tan 2 p g VTW J ' 
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where TVj 3 , mu 3 are the third generation soft mass parameters and A represents the scale at 
which supersymmetry breaking is mediated. For simplicity, we assume the decoupling limit 
such that rrih denotes the physical Higgs boson mass and the mass scales in the logarithmns 
have been replaced by the TeV scale. We will also neglect A-terms (i.e. At), which are 
generally predicted to be small for gauge mediated supersymmetry breaking. Even if they 
are not (as occurs in some models [8,21-23]) our conclusions are not affected as all bounds 
apply to the lightest, physical stop mass. 

Requiring that the fine tuning is no worse than ten percent (i.e. A < 10) and assuming 
that A is not much bigger than 10-100 TeV leads to the following constraints on the soft mass 
parameters: 

mQ 3 ,ma 3 < 500 - 700 GeV, M 3 < 1100 - 2200 GeV, (5) 

fi < 300 GeV, M 2 < 1500 - 2000 GeV. (6) 

All sfermions other than the stops and the left handed sbottom can be decoupled since they 
do not significantly affect the fine tuning in the Higgs sector^] Hence a common scale will be 
assumed for the remaining soft mass parameters 

m Q 1 y m «l,2> m <f li2 , 3 ' m Ll,2,3' m ei,2,3> m i%2,3 = M . ( 7 ) 

where M > few TeV. This leaves the minimal sfermion spectrum required for naturalness. 
In the gaugino sector we first assume the standard one loop relationship 

3 Mi = M 1= M 1 

5 OL\ Ct2 Q!3 ' 

i.e. gaugino unification, continues to hold. Combined with the existing gluino mass bound 
from naturalness this implies 

M i = \^r£r\ M * £ 160 " 320 GeV ( 9 ) 

3 a 3 (M z ) 

which will be important in determining the identity of the NLSP. In the Higgs sector we fix 
rrih = 125 GeV. A number of examples with an NGM spectrum [T2"|ll4pi5pi7pi8] already con- 
tain a mechanism to raise the Higgs mass. However, it should be stressed that any mechanism 
which raises the Higgs mass may also contribute to the amount of fine tuning. 

Little parameter space remains for stop or sbottom NLSPs. Searches for direct stop 
production [2"7W3"T] exclude stop masses between 220 and 500 GeV, and between 110 and 165 
GeV. In addition, left handed stop or sbottom NLSPs below 350 GeV and a right handed 
stop NLSP below about 200 GeV were excluded in Ref. [5] (although a stop NLSP close to 
the top mass is not yet excluded). Since eq. jSJ) implies that M\ < M 2 < M3 there are two 
further possibilities. The first is a Higgsino-rich chargino, requiring \fj>\ < \M\\ < 160 — 320 
GeV. This scenario turns out to be somewhat non-generic [32] so we do not focus on it here 
(although chargino NLSPs outside of an NGM context have been investigated in Ref. [33]). 
By far the most common NLSP is a bino or Higgsino-rich neutralino, also lighter than about 



1 Sleptons only have a small effect on LHC phenomenology hence our results also apply to models where 
they are not decoupled. The exception is when the NLSP is a slepton, typically a stau, which is common if 
the entire third generation is kept light. Generalising stau searches [24H26] to include light stops is therefore 
a well motivated extension to this work. 
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Figure 1: A typical spectrum for natural gauge mediation. 



160 — 320 GeV. NGM models with a Higgsino-rich neutralino NLSP decaying to Z bosons 
were analysed in Ref. [33], so we focus on bino-rich neutralinos, which decay to photons and 
gravitinos. A typical spectrum is given in Figure [1] 

Relaxing the constraint imposed by eq. ([5]) the spectrum may remain similar or may be 
quite different. For example, it is well known that many explicit models only generate gaugino 
masses at two loops. This further suppresses M\ so the NLSP is still a bino-rich neutralino. 
On the other hand, one can easily construct more involved models where the gaugino mass 
hierarchy implied by eq. ([8]) is completely disrupted, whereupon alternative NLSPs become 
possible. To constrain such models one would need to repeat the relavant analysis in Ref. [6], 
but with a light third generation. Generically one would expect to find stronger bounds on 
gluino and NLSP masses, but weaker bounds for the stop mass relative to those found for 
degenerate squarks. 

In order to constrain the NGM framework we take a purely phenomenological approach 
using simplified models. This enables robust bounds to be placed on physical masses inde- 
pendently of model details. Since fine tuning requirements force all NGM models to contain 
light Higgsinos, light stops and light gluinos, one expects both coloured and electroweak spar- 
ticle production. We consider separate simplified models for each process. In practise both 
processes contribute in any NGM model, hence our bounds are always conservative. 

For coloured production we decouple all sparticles other than the gluino, the right-handed 
stop and the bino. Hence sparticle creation proceeds through gluino or stop pair production. 
Including other light sparticles would strengthen our final bounds. The bino must be light 
enough such that the NLSP is mostly bino (so decays to photons) but otherwise its mass has 
only a weak kinematic effect on the signal strength. This can be seen in Refs. [35J136] and we 
have verified that the situation is the same here. We thus fix the bino mass at 100 GeV in 
this scenario. For electroweak production we decouple all coloured sparticles and keep only 
the Higgsinos and bino in the electroweak sector. Sparticle creation now proceeds through 
neutralino/chargino pair production. Naturalness does not forbid heavy winos and including 
them among the light states again leads to stronger bounds. 

3 Exclusion limits 

Neutralino NLSPs decaying to photons lead to a very distinctive signature: two high energy 
photons and missing transverse momentum from the gravitinos. This channel has a low back- 
ground and is the subject of recent updates from the ATLAS [35] and CMS [36] experiments. 
Only the ATLAS diphoton search is considered here, on the basis that CMS and ATLAS have 
a similar reach. We also considered the earlier, 1 fb _1 ATLAS search [37], which has softer 
kinematic requirements, but found that it offered no additional constraints on NGM. 

Three signal regions were defined in the ATLAS diphoton analysis. Of these the first two 
are the most constraining for coloured production and the third for electroweak production. 
At least two isolated photons with px > 50 GeV are required, for which the energy in a cone 
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Region 


Rl 


R2 


R3 


Photons (E T > 50 GeV) 


> 2 


> 2 


> 2 


Ef iss (GeV) 


> 200 


> 100 


> 125 


A0( 7 ,P? iSS )min 


> 0.5 




> 0.5 


H T (GeV) 


> 600 


> 1100 




Expected background 


0.10 ± 0.03 ± 0.07 


0.36 ± 0.05 ± 0.27 


2.11 ± 0.37 ± 0.77 


Observed events 
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95% CL upper limit 


3.1 


3.1 


4.9 



Table 1: Selection cuts for the three ATLAS diphoton search signal regions. A<^(7,p™ lss ) m j n 
is the smallest of the azimuthal separations between the missing momentum p^ lss and the 
momenta of the two leading photons in the event. The total visible transverse energy Ht is 
the scalar sum of the transverse momenta of the jets, leptons and two leading photons in the 
event. 

of AR = 0.2 surrounding the photon's deposition in the calorimeter must be less than 5 GeV. 
Table [1] summarises the further selection cuts for each signal region, the results of the ATLAS 
search, and the ATLAS derived limits on the number of events attributable to new physics. 
These results may be used to constrain the NGM parameter space if we evaluate the signal 
expectation for candidate NGM models. 

We use SUSY-HIT 1 . 3 [38] to produce mass and decay spectra for a given set of model 
parameters, Pythia 6.24.26 [3U] to generate 10,000 Monte Carlo events for each point, and 
a custom version of Delphes 1 . 9 [30] (with photon isolation added) to provide a fast ATLAS 
detector simulation. The total supersymmetric production cross section is calculated at next- 
to-leading order using PR0SPIN0 2.1 [IT]. A cone-based overlap removal procedure is used 
to avoid double counting particles that are reconstructed as more than one object (e.g. an 
electron and a jet). 

To approximate the ATLAS limit setting procedure we follow the approach used in 
Ref. [32]. We use the published ATLAS limit in the m^o — m g plane to calibrate system- 
atic error parameters in a simplified model of the ATLAS likelihood function, then use this 
likelihood function to generate limits in the NGM parameter space. We have confirmed that 
our procedure reproduces the ATLAS 95% confidence level exclusion contours and, where 
modest discrepancies are encountered, we tune our parameters to ensure that our results are 
more conservative. See Figure [2] for details. 

In our exploration of the NGM parameter space, we fix tan f3 = 2, M = 2.5 TeV and 
cr^o < 0.1 mm (to ensure prompt neutralino decays) throughout. For coloured production we 

fix ttIq 3 , fi and M2 at the high scale M, M\ = 100 GeV, then scan over m« 3 (approximately 
equal to when nin 3 > nit) and M3 (i.e. the gluino mass nig). For electroweak production 

we instead fix m^, niu 3 , Mi and M3 at the high scale M, then scan over Mi (approximately 
equal to m^o when Mi < fi) and \i (approximately equal to m-± when M\ < fj,). 

The resulting limits are shown in Figure [3l For coloured production the gluino mass limit 
approaches ~ 1100 GeV for heavy stops, in agreement with the published limits from ATLAS 
and CMS [331136] (note that varying the bino mass does not significantly change the limits as 
can be seen from Figure [2]). Consequently one must accept a fine tuning of at least 10 — 40% 
due to two loop Higgs mass corrections from gluinos. Stop masses less than ~ 580 GeV 
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Figure 2: Comparison between the Delphes and ATLAS 95% exclusion limits in the m^o — m g 
plane. The Delphes limit is obtained by taking the union of the Delphes limits for each 
signal region, with the limits tuned to match the ATLAS results. 



are excluded (the effect of including weak production as well would raise this limit further), 
corresponding to a fine tuning of at least 8 — 17% due to one loop Higgs mass corrections 
from the top Yukawa coupling. These limits do not degrade at low stop mass as photons from 
neutralino decays typically remain hard enough to pass all cuts for m^o = 100 GeV. 

Limits on weak production reveal that a significant slice of the m^o — m~± plane is also 

Al Xl 

excluded. The limit does not extend to neutralino masses below 50 GeV as too few events in 
these models pass the missing energy requirement of the ATLAS search. Nor does it extend 
to m^o ~ m~± where the branching ratio of neutralinos to photons drops off. One could in 

Xi Xi 

principle have models with m^o < 50 GeV and a small value of [i (and hence low fine tuning) 
although, if one insists on gaugino unification, this is difficult to achieve due to limits on the 
gluino mass from coloured production. If m { o > 50 GeV one must have m-± > 300 — 470 
GeV and a fine tuning of at least 4 — 9% from the ^i-term. 

Combining both limits one can search for the most natural NGM model (with a neutralino 
NLSP) not yet excluded by the LHC. If gaugino unification is assumed both charginos, from 
the Higgsino and wino, are light. Constraints from electroweak production are therefore 
severe. Choosing 

m-g > 2200 GeV, u. > 1000 GeV, (10) 

avoids both limits and yields a fine tuning of at least 1% from the /x-term. The fine tuning from 
the top Yukawa coupling is subdominant. The constraints are relaxed if gaugino unification 
is not assumed such that the wino is allowed to be heavy and the bino light. Then 

m g > 1600 GeV, n > 400 GeV, m^, > 800 GeV, (11) 

avoids both limits with a fine tuning of at least 5% from the //-term and 5 — 9% from the 
top Yukawa coupling. As expected, the fine tuning in a realistic spectrum is worse than that 
suggested by the simplified models alone. 
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Figure 3: The 95% confidence level exclusion contours. Left: Limits on coloured production in 
the mjj — mg plane (note that the limit does not degrade at low stop mass for m^o = 100 GeV). 
Right: Limits on electroweak production in the m^o —m~± plane. Both production processes 
are active in any given model so all limits are conservative. 



Finally, we note that the current ATLAS search has not been optimised for direct stop 
production. A more dedicated search using 6-tagging and reoptimised selections on E™ ms and 
Ht could increase the reach, and we strongly encourage effort in this area. 

4 Conclusion 

In summary, using a 4.8 fb" 1 ATLAS diphoton search, we have placed mass limits on simplified 
models of natural gauge mediation with a neutralino NLSP for coloured and electroweak 
sparticle production. Because both production mechanisms are active in actual realisations 
of natural gauge mediation, our bounds are conservative. Top squarks with a mass below 
~ 580 GeV are excluded, as are gluinos with a mass below ~ 1100 GeV. Assuming the 
neutralino is heavier than 50 GeV means that charginos lighter than ~ 300 — 470 GeV are 
ruled out. Otherwise a neutralino lighter than ~ 50 GeV may allow our bounds to be evaded. 
This places the first nontrivial constraints that tests the naturalness of this class of models. 
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